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^ . Abstract 

The variational approach to the Hamihon formulation of Yang-Mills theory in Coulomb 
^ ' gauge developed by the present authors previously is applied to Yang-Mills theory in 2-|- 1 

CN , dimensions and is confronted with the existing lattice data. We show that the resulting 

Dyson-Schwinger equations (DSE) yield consistent solutions in 2 -|- 1 dimensions only for 
(S| ' infrared divergent ghost form factor and gluon energy. The obtained numerical solutions 

of the DSE reproduce the analytic infrared results and are in satisfactory agreement with 
the existing lattice date in the whole momentum range. 
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^ : 1 Introduction 

Recently, there has been a renewed interest in the formulation of Yang-Mills theory in Coulomb 
gauge, both in the continuum theory [H [2], [3l H] and on the lattice [5], El [3, El E]. In the 
continuum formulation the Hamilton approach turns out to be very appealing, in particular, 
for the description of the confining properties of the theory |3J. The reason is that in Coulomb 
gauge Gauss' law can be explicitly resolved resulting in a static potential between color charges 
|10j . Recently, several papers have been devoted to a variational solution of the Yang-Mills 
Schrodinger equation in Coulomb gauge [Zl [3l H]. In particular, in refs. [Sj H] the present 
authors have developed a variational approach to Yang-Mills theory in Coulomb gauge which 
properly includes the curvature of the space of (transversal) gauge orbits. Using a Gaussian 
type of ansatz for the Yang-Mills vacuum wave functional minimization of the energy density 
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results in a set of coupled Dyson- Schwinger equations which can be solved analytically in the 
infrared [H [TT] and has been solved numerically in the full momentum range [51 EI] in Z) = 3 + 1 
dimensions. An infrared divergent gluon energy and a linearly rising static quark potential was 
found [3l [21], both signaling confinement. Furthermore, within this approach the spatial 't 
Hooft loop [22], [23] was calculated in the vacuum state and a perimeter law was found [20] , 
which is the behavior expected for this disorder parameter in a confinement phase. In the 
present paper we apply this approach to 2 + 1 dimensions and confront the results with the 
existing lattice data p5]. 

The 2 + 1 dimensional Yang-Mills theory is interesting in several aspects: When a Higgs field is 
included, it represents the high temperature limit of the 3 + 1 dimensional Yang-Mills theory, 
thereby the temporal component of the gauge field Aq becomes the Higgs field in the dimensional 
reduced theory. In many respects the 2 + 1 dimensional theory is easier to treat than its 3 + 1 
dimensional counter part, in particular, much larger lattices can be afforded in 2 + 1 dimensions. 
This will be crucial for a comparison of the continuum results with the lattice data, since the 
lattices affordable in 3 + 1 dimensions are by far too small to allow for a reliable extraction of the 
infrared properties of the Greens functions [121 [13]. addition, 2 + 1 dimensional Yang-Mills 
theory is super-renormalizable. 

The balance of the paper is as follows: In section 2 we briefly summarize the essential ingredients 
of the Hamilton approach to Yang-Mills theory in Coulomb gauge and of the variational 
solution of the corresponding Yang- Mills Schrodinger equation [3]. In section 3 the Dyson- 
Schwinger equations resulting from minimizing the energy density are studied in the ultraviolet 
and their renormalization is carried out. Some exact statements on their solutions are given 
in section 4, where we also solve these equations analytically in the infrared. In section 5 
we present our numerical results and compare them with the existing lattice data. A short 
summary and some concluding remarks are given in section 6. 



2 The Hamilton approach to Yang-Mills theory in Coulomb 



Canonical quantization of Yang-Mills theory is usually performed in Weyl gauge = to avoid 
the problem arising from the vanishing of the canonical momentum conjugate to Aq} The prize 
one pays by choosing Weyl gauge is that one looses Gauss' law as equation of motion. To ensure 
gauge invariance one has to impose Gauss' law as a constraint on the wave functional. Instead 
of using gauge invariant wave functionals [HI [151 [iH] , it is simpler to explicitly resolve Gauss' 
law by fixing the residual (time-independent) gauge invariance, and Coulomb gauge diAi = is 
a particularly convenient gauge for this purpose. After resolving Gauss' law in Coulomb gauge 
the Yang-Mills Hamiltonian reads [TO] 



gauge 



H 




(1) 



^An alternative to the Weyl gauge is the light-cone gauge Aq ± Ai = 0. 
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where 

is the canonical momentum operator conjugate to the transverse gauge field A-^. Furthermore 

J[A^] = det (^-dA{A^)) (3) 

is the Faddeev-Popov determinant with D{A) = d + gA , A = A'^Ta , {^a^ — being the 
covariant derivative in the adjoint representation. Furthermore 

B = ^e''[D,, D^] , A = 5. + igAi (4) 

is the magnetic field, which is a scalar in 2 + 1 dimensions (ejj = — e^j, ei2 = 1). Finally 

F»''(x,x') = (xa|(-A9.)-^(-92)(-4-9,)-i|x'6) (5) 

is the non-Abelian Coulomb propagator which mediates a static interaction between the color 
charge density of the gluons 

p'^(x) = -i,^"^(x)n,^^(x). (6) 

In the presence of external color charges, for example in the presence of quarks, their charge 
density has to be added to the gluon charge density. The kinetic term in the Hamilton has 
the form of a (variational extension of the) Laplace-Beltrami operator in curved space with the 
Faddeev-Popov determinant ([3]) corresponding to the determinant of the metric of the space of 
transversal gauge orbits. The Jacobian ([3]) also enters the scalar product in the Hilbert space 
of the Yang-Mills wave functionals in Coulomb gauge 



(^i|0|^2) = j VA^J[A^]<^l[A^]0[A^,Ii^]^2[A^] . (7) 

We will solve the Yang-Mills Schrodinger equation by the variational principle 

E = min (8) 



using the following ansatz for the vacuum wave functional [3], [1] 

^A^] = =Arj^[A^]-°exp J d?x j rfV4^"(x)cu(x,x04^"(x') 



(9) 



where uj{x, x') is the variational kernel, which by translational and rotational invariance of the 
vacuum depends only on Ice — x'\, and, by the isotropy of color space, is independent of color. 
The ansatz ([9]) with a = | is motivated by the wave functional of a particle moving in a s-state 
in a spherically symmetric potential. In principle, a could be used as a variational parameter 
to minimize the energy density. However, it turns out that up to two loops in the energy. 
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stationarity of the energy density with respect to u, i.e. SE/Su = 0, imphes also stationarity 
with respect to a, i.e. dE/da = [4]. Thus, we are free to choose a for convenience, and 
as in ref. |3j we will choose a = |, which removes the Faddeev- Popov determinant from the 
integration measure in eq. ([7j). Furthermore with the choice « = | the gluon propagator is 
given by 

{^Ai{x)Af{ym = ^U,{x)u-\x,y) , (10) 
so that the Fourier transform u){k) has the meaning of the gluon energy. 

The calculation of the vacuum expectation value of the Coulomb Hamiltonian in the state ([9]) 
proceeds in the same way as in 3+1-dimensions and we just quote the result. We find for the 
kinetic energy 

4 J uj[k) 



the potential energy 

E, = ^^^1^5(2)(o) j d'k 



^ ^c{Nl-i) r dW 1 / _ (kkr . 

+ 16 3d ^^)J (27r)2 ^(kMkO V k%'2 ' ^ ' 



and for the Coulomb energy 

^ iVc(iVg-l) ,(2).p^ f d'kd'k' {Wfd{k-krf{k-k') 
16 ^'J (27r)2 k^k'' {k-k'Y 

{[uik) - xm - Hk') - x{k')] f 



uj{k)uj{k') 



(13) 



where 5^(0) = \^/(27r)^ with V = J d'^x being the 2-dimensional spatial volume. Furthermore 
the quantity 

= -Nf^i {jaP^MfW)), ^^^^ 

is referred to as "curvature" of the space of gauge orbits. Introducing the ghost propagator by 

-1 - ^id{-d^ 

9 

with d{k) being the ghost form factor, the curvature ( IT^ can be expressed as 

Nc f d\ ( (kq)2\ rf(k-qKq) 



G = {^\(-d,DA \^) = (-d'y^^ ^, (15) 
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Finally, f{k) denotes the Coulomb form factor, which is defined by 



{^\ (^-dAy\-d') (^-dAy'm = G{-d')fG. (17) 

Minimization of the energy density {E^ + Ep + Ec)/V results in the so-called gap equation, 
which can be cast into the form of a dispersion relation of a relativistic particle (gluon) 

a;(k)2 = k2 + x(k)^ + /.(k) + /°. (18) 

Here I'^ is an irrelevant constant which arises from the gluon tadpole. Furthermore, the quantity 
Iuj{k), which arises from the expectation value of the Coulomb term ( fT3l) . can be expressed as 



/ (k): I' "^'^ (kq)^rf(k-q)V(k-q 



2 J (27r)2 k2q2 (k - q)^ 

[^(q)-x(q) + x(k)]'-c.(k)^ 



(19) 



In principle G or d , x ^^^^ f ^ire defined uniquely once the trial wave functional "^{A) is 
fixed. However, the exact evaluation of these expectation values (even with the above chosen 
trial wave function) is not feasible and one has to resort to further approximations. We shall 
adopt here the same approximation as in ref. [3] used in 3 + 1 dimensions, approximating 
the full ghost-gluon vertex by its bare one. This approximation has been justified in Landau 
gauge [17] and, in fact, has received recently strong support by lattice calculations [18j. In 
this approximation the ghost and Coulomb form factors satisfy the following Dyson-Schwinger 
equations 

1 - /d(k) , (20) 



d{k) 

W = / (1^ (l - ^) (k-''q)2^(q)' ^2^) 
/(k) = l + /;(k), (22) 



r -^O f A _ (kq)2\ rf(k-q)2/(k-q) 

As described in ref. [3] to the considered order (bare ghost-gluon vertex and 2-loop order in 
energy) the Coulomb form factor can be put f{k) = 1. 

In Z) = 2 + 1 the coupling constant g has the dimension of the square root the momentum and 
it is convenient to use the coupling constant to rescale all quantities by suitable powers of g to 
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render them dimensionless. The couphng constant then disappears from the Dyson-Schwinger 
equations. Denoting the dimensionless quantities by a bar we have 



T k _ - uj{g^k) _ - xig^k) 
= - , uik) = , xik) = 

m = frk) = fi9'k). (24) 

9 

In the following we will skip the bar and, unless stated otherwise, all quantities will be under- 
stood as the dimensionless ones. 

In section 14.11 we prove that for any solution of the coupled Dyson-Schwinger equations (fT8l 
[2U]) in D = 2 + 1 the ghost form factor and the gluon energy are infrared divergent 

d-\k = 0) = (25) 
uj~\k ^0) = . (26) 

The first relation is the so-called horizon condition which had to be imposed ad hoc in D = 3 + 1 
but is a strict consequence of the Dyson-Schwinger equations in D = 2+1. The second condition 
fl2B]) signals gluon confinement. In addition, we prove in appendix A that when the curvature 
x{k) is omitted, the coupled Dyson-Schwinger equations in D = 2 + 1 do not allow for a 
consistent solution. This is again different from D = 3 + 1 where the DSE allow for solutions 
with infrared finite uj{k) when the curvature x{k) is ignored. 



3 Ultraviolet behavior and renormalization 

In the following we will investigate the ultraviolet behavior of the solutions of the Dyson- 
Schwinger equations and perform their renormalization, thereby following the procedure pre- 
sented in refs. [3], [201 ^^r the D = 3 + 1-dimensional case. Since most of the considerations 
will parallel the 3 + 1-dimensional case we will be very brief. 
For the uv-analysis it is sufficient to use the angular approximation 

h{\k - q\) = h{k)Q{k -q) + h{q)e{q - k) . (27) 

The angular integral in the ghost Dyson-Schwinger equation and also in the curvature then 
becomes trivial J dip sin^ ip = t[ and we obtain for the corresponding momentum integrals 



h{k) 



Nc 
Sir 



X{k) = 



Nc 
Sir 



d{k) 



d{k) 



dq 



uj{q) 



+ / dq 



dqqd{q) + dq 



d{q) 
quj{q) 

d{q? 



(28) 



(29) 
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The asymptotic analysis of the Dyson- Schwinger equations is simphfied by taken the derivative 
with respect to the external momentum 



Ncl_ 

Ncl_ 
Stt k^ 



d'{k) - 2 



d'{k) - 2 



d{k) 



d{k) 



dq 



uj{q) 



dqqd{q) . 



(30) 



(31) 



As we will see below, the remaining momentum integrals are ultraviolet convergent^. The 
derivative of the ghost Dyson-Schwinger equation ( |20l) yields 



d'{k) 



NcR{k) 



_d{ky Stt F J An k F 



Nc d{k) R{k) 



where 



K 

R{k) = j dq 



uj{q) 



and the derivative of the curvature (fT6|) yields 



x'ik) 



Ncl_ 
8tt A;2 



d'{k) - 2 



d{k) 



Sik) 



where 



S{k) = J dqqd{q) . 



We now discuss the ultraviolet behavior of the relevant quantities. 



(32) 



(33) 



(34) 



(35) 



3.1 Ultraviolet behavior 

For large momenta the k"^ term on the r.h.s. of the gap equation (|T8|) dominates [3] and the 
gluon energy uj{k) behaves like 

ijj{k) —>■ , k ^ oo (36) 

in accordance with asymptotic freedom. This behavior will in fact be confirmed by the numerical 
solutions presented in section 6. Assuming (!36|) . we will investigate the uv-behavior of the 
remaining quantities: the ghost form factor d{k), the curvature x{k) and the Coulomb form 
factor f{k). 

^The ghost integral Id(k) is uv-finite as will be shown later. 
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For large k —>■ oo the fc-dependence of the integral (l33l) is independent of the infrared behavior 
of uj{k) and with eq. fl5^ we find 



R{k) = k , k ^ oo 
and the derivative of the ghost DSE ( l32l) reduces to 



d'{k) 
d{ky 



Ncdjkf 
Stt k 



Ncdjk) 
4:Tc k'^ 



(37) 



(3^ 



To solve this equation, let us assume for the moment ^^^^ -C 1 for k oo. Then the differential 
equation (138|1 reduces to 



d'{k) Nc d{k) 



whose solution is given by 



d{ky 



d{k) 



4:71 



L _ 1 

-2 Stt k 



(39) 



(40) 



where c is an integration constant. Indeed, for large k ^ oo this solution satisfies <^ 1. 
To determine the integration constant c, we consider the asymptotic behavior of the integral 
Id{k) for k oo, where uj{k) ~ k and d{k) ~ c. This yields 



c dq 



quj{q) 



Ncc 
Ait k 



(41) 



Since Id{k oo) —>■ 0, we obtain from the DSE of the ghost form factor d{k oo) —>■ 1, which 
fixes the integration constant in eq. fHDj) to c = 1, so that the asymptotic form of the ghost 
form factor becomes 



d{k) 



1 



1 _ 1 

Stt fc 



— >■ OO . 



(42) 



Accordingly we find for the unsealed form factor'^ d{k) = gd i-^j the asymptotic form 



d{k) 

From the Swift relation [19], [3] 



1 _ ilC 31 



k ^ OO . 



(43) 



1 d 1 



(44) 



^ Below we denote the dimensionful (unsealed) quantities by a "tilde". 
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we find for the Coulomb form factor 

f{k) = ^ ^ = d{k) , k^oo, (45) 
1 _ 

87r fc 

which obviously has the same asymptotic form as the ghost form factor. 

With the asymptotic behavior of the ghost ( l42l) . we find for the derivative of the curvature ( |34l) 

X'{k) = -^l,S{k) ,k-.oo (46) 



S{k-.oo)^^, (47) 



x\k) = ^l,k^oo, (48) 



X{k) ~ In ( ^ ) , A; ^ oo . (49) 



X{k) li4^y_2^o (50) 



and for the integral S{k) ([35 



so that 



I.e. 



Accordingly, the ratio 



uj{k) k \fi 

vanishes in the ultraviolet implying that the space of gauge orbits becomes asymptotically fiat 
in accordance with asymptotic freedom. 

3.2 Renormalization 

With the above obtained uv-behavior (see eqs. (136|) . (H2l) and f H5|) ) the integrals Id{k) fl2T]) and 
If{k) (123|) are uv-convergent. Thus, contrary to the D = 3 + 1 dimensional case in D = 2 + 1 
the DSE for the ghost and Coulomb form factors do not need renormalization. What needs, 
however, renormalization is the curvature x{k) and the gap equation. The renormalization of 
these two equations is carried out in exactly the same way as described in refs. [20] and |3], 
i.e. basically by subtracting this equations at a renormalization point fi, which leads to the 
renormalized equations ^0] 



X{k) = x{l^) + m,m = Ixik)-W (51) 

iv\k)-x\k) = A;2 + eo + A/(f)(A;) + 2x(A;)(e + Aii,i)(A;)-Aii^)(0)) . (52) 

Here xif^) and 

e = x(/i) + ii^^(0) (53) 
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are finite renormalization constants. Furthermore, we have introduced the abbreviation 



^il^\k) = i^:;\k) - (54) 



with 



Ti^)ik A^-^^ r (i.^^^ rf(k-q)V(k-q) [^(q) - x(q)]" - k(k) - x(k)]" 



(55) 



Note, since the ghost equation needs no renormahzation, there are only three renormahzation 
constants ^o, i and x(/^) and furthermore, the solutions of the coupled ghost and gluon DSEs 
do not depend on (Note the integrals ( l55i) depend only on the finite quantity x(^) but 

not on the renormalization constant x{^^) [20] •) One of the two independent renormalization 
constants, ^o, is used to implement the horizon condition d~^{k — 0) = 0. As will be shown in 
the next section any consistent solution of the coupled DSEs does satisfy this condition. The 
remaining renormalization constant C, can be chosen at will and determines the infrared limit 
of the wave functional as will be shown at the end of the next section. 



4 Analytic results 

The Dyson-Schwinger equations arising from the variational solution of the Schrodinger equa- 
tion have in principle the same form inD = 2 + lasinD = 3 + l dimensions. However, due 
to the fact that the 2 + 1 dimensional theory is superrenormalizable some rigorous properties 
of the solution of the Schwinger- Dyson equations can be derived which are not accessible in 
3 + 1 dimensions. Below we shall derive some rigorous properties of the solutions of the DSE 
and discuss their physical implications. 

4.1 General results 

Consider the ghost DSE (1201) . Since the integral Id{k) (12 ip is convergent in = 2 + 1 (unlike 
in D = 3 + 1), the ghost DSE needs no renormalization and thus no renormalization constant 
is introduced by this equation. It is then not difficult to prove the following statement: 

If d{k) is a continuous function in A; G [0, oo), it satisfies 

d{k) > 1 . (56) 



We prove the statement by reductio ad absurdum: Assume d{k) < for all k G [0, oo). Then, 
since u}{k) > by normalizability of the wave functional it follows from (12T]) Id{k) < and thus 
from the ghost DSE ([20]) 

d{k) > (57) 



h{k) 



*Some observables like the 't Hooft loop (in 13 = 3 + 1) do, however, depend on yf/i) |20j. 
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in contradiction to the assumption. Hence, d{k) cannot be negative everywhere. Assume now 
there exist some momentum k' for which d{k') < 0. Since d{k) as a solution of the DSE, can be 
assumed to be continuous and as shown above is not everywhere negative, it must have at least 
one zero, say ko, i.e. d^k^) = 0, where d{k) changes sign. By the ghost DSE, eq. (120|) Id{k) 
has to be singular at k ^ ko and change sign at k = ko, too. However, changing integration 
variable in the (convergent!) integral Id{k) (!2T|) we find 



Nc 

Id {ko ± eko) = — 



_fq_ 



1 



(fcp ■ {q ± eko)y 
kl {q ± ekof 



d{ko 



{ko -q) uj{q- eko) 



{51 



and by the positivity of uj{k), Id {ko + eko) and Id {ko — eko) have the same sign in contradiction 
to the above assumption. Thus, d{k) > holds for all k. Then from eq. fl^Tl) follows Id{k) > 
and by the ghost DSE fl^Ul) d{k) = (1 — Id{k))~^ > 1. Note, that (jSSD is also in agreement with 
the asymptotic uv-behavior found in the previous section. 

Assume now that the ghost form factor is bounded from above, i.e. there exist some upperbound 
M > 1 such that d{k) < M for V/c G [0, oo). According to (156!) d{k) is then restricted to the 
intervall 1 < d{k) < M and the integrand in x(/c) (fT6l) is positive definite. Therefore, replacing 
d{k) in the curvature (ITB]) by its upper and lower bound, we obtain an upper and lower bound 
to x{k) 



MH > x{k) >I,I 



Nc 
2 



d\ 
{2ny 



1 - 



k2q2 



q) 



(59) 



where we have introduced a momentum cutoff A since the integral is uv-divergent. The angular 
integral can be done analytically^ 



TT 

/ 



dip 



sin^ if 



/j2 _|_ g2 _ 2kq cos 



TT 

2 



1 
1 



q < k 
k < q 



(60) 



yielding 



N, 



c 



87r2 



7 



dqq 



TT 



Nc 
Svr 



1 , A 
- + ln- 

2 k 



(61) 



Inserting this result into eq. fl59|l . we find for x{k) = x{k) — x(/^) the bounds 



-^MHn^>x{k)> 
Sir jjL 



Svr jjL 



(62) 



Note that this relation holds for all k. In particular, this relation shows that x(/c) is infrared 
divergent 



m 



In — , /c 



. 



(63) 



^Note that this exact relation is also reproduced by the angular approximation. 
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For an infrared divergent x(^) the gap equation fl5T]) reduces in the infrared to (see also ref. 



liniu;(A;) - x{k) = ^ 



fc-»0 



(64) 



where ^ is the renormahzation constant introduced in eq. fl53l) . Eq. fl64|) shows that uj{k) has 
the same infrared behavior as x{k)- Thus, if the ghost form factor is bounded 1 < d{k) < M 
the gluon energy is logarithmically infrared divergent u(k) ~ — In-. It is now not difficult to 
show that with such an infrared behavior of uj{k) the ghost Dyson-Schwinger equation does not 
possess a solution. To show this let us assume that there exist an e > such that 



u{k) < a ( - In^ ) , VA; G (0,e) 



(65) 



with some positive constant a, which includes the case 
Dyson-Schwinger equation (^11) 



Consider the integral in the ghost 



Id{k) 



27r 



C 



87r2 



dq 



uj{q) 



dip sin if 



2 . c?(k - q) 





Since d{k) > 1, see eq. fISBl) . we obtain the following estimate 

oo 27r e 



(k-q) 



(66) 



27r 



Id{k) > 



N, 



c 



dq 



uj{q) 



dip sin^ Lp ■ 



(k-q) 



> 



c 



dq 



uj{q) 



dp) sin^ p> 



(k-q)2- 



(67) 



The angular integral can be done exactly using eq. fl60l) resulting in 

k e 



Id{k) > 



c 



dq 



q n 
uj{q) fc2 



+ / dq 



q TT 
uj{q) g2 



c 



>^ dq 



quj{q)' 



Inserting here eq. (|65|1 we obtain 



N, 



h{k) >^-^ I dq 



N, 



c 



ga(— In^) Sira 



In 



In 



In^ 



oo . 



Thus we find 



Id{k) ^ oo. 



(68) 



(69) 



(70) 



and from the ghost Dyson-Schwinger equation ( !20|) follows d{k 0) = 0, which is in con- 
tradiction to the rigorous result ( l56l) . Thus we have shown that the coupled ghost and gluon 
Dyson-Schwinger equations do not allow for a ghost form factor which is bounded from above. 
We now show that ghost form factor d{k) is a monotonously decreasing function of k, i.e. 



d'{k) < 



(71) 
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for all finite k. 

As shown above in the ultraviolet analysis d'{k) < for k — > oo. Assume now, as we lower k, 
at some finite k = ko, d'{k = ko) = 0. Then from the eq. fl5Ul) follows I'Jyk^) < 0, which is in 
contradiction to the ghost DSE from which follows 

Thus flTTl) holds in the whole momentum range. Since fl56|) d{k) > 1 and d{k) is monotonously 
decreasing in the whole momentum range k > and furthermore d{k) must not be bounded 
from above, it follows that d{k) is infrared divergent, i.e. 

d-\k = 0) = , (73) 

which is the horizon condition. 

This is different from the 3 + 1 dimensional case where solutions to the Dyson-Schwinger 
equation exist with an infrared finite ghost form factor. 

4.2 Infrared analysis 

The DSE can be solved analytically in the infrared completely analogous to the D = 3 + 1 
dimensional case. For this purpose we make the following power ansatze in the infrared 

-(^) = ^ ' dik) = ^ , x{k) = ^ . (74) 

We will first resort to the angular approximation used already in the uv-analysis. Later on we 
will present the results obtained without resorting to the angular approximation. 
With the infrared ansatze (ITill . we find for the integrals defined by eq. (133|1 and (l35l) 

R{k 0) = ^j^k'^'^ , S{k ^ 0) = j^k'-^ . (75) 
Al + a 1 ~ P 



From the derivative of the ghost DSE (!72|) we obtain 

A Nc (3 + 2 



Stt p{a + 2) 

which implies 



, (76) 



a = 2(3 (77) 



and 

A Nc (3 + 2 



(7J 



52 Stt 2(3{(3 + 1) ' 

In an analogous fashion we obtain from the derivative of the curvature x'{k) = I^{k) 

^ = ^ ^ + ^ F-^/^ (79) 
52 Stt 7(2-/5) ^ ^ 
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resulting in 



and 



From eq. (1771) and (IHO!) follows 



7 = 2/3 (80) 

C _Nc P + 2 

Stt 2/3(2 -(3)' ^ ' 

a = 7 = 2/3 , (82) 



showing that uj{k) and x(^) have the same infrared exponents just like in the D = 3 + 1 
dimensional case (ref. [3]). In fact, a = 7 follows already from the infrared limit of the gap 
equation fl64p . From this equation in addition follows that not only the infrared exponents but 
also the prefactors of both quantities have to coincide, i.e. 

A = C . (83) 

Dividing eq. flTHj) by flHTl) and using fl82|) we obtain 

A (2-/3) 
C (1 + /3) 

and A = C 



(84) 



/? = ^ , a = 7 = 1 . (85) 

As we will see in sect. 6, this infrared behavior yields a linearly rising static color Coulomb 
potential provided we approximate the Coulomb form factor f{k) by its leading term f{k 
0) = 1 which is correct to the order considered in the present paper. 

The infrared analysis of the Dyson-Schwinger equations can be also carried out without resorting 
to the angular approximation. In fact, in ref. [IT] the infrared analysis was carried out for 
arbitrary dimensions. In that case one finds from the ghost DSE the following sum rule for the 
infrared exponents 

a = 2(3 + d-2 (86) 

due to the non-renormalization of the ghost-gluon vertex. As shown in ref. [TT] this sum rule 
guarantees that x{^) ^iid uj{k) have the same infrared exponent a = 7 as already found above 
in the angular approximation, and in agreement with the infrared limit of the gap equation 
(16^ . Eq. (IMIl together with the ghost DSE in the infrared limit can be solved analytically for 
the infrared exponents yielding [TT] 

13 = 0.4 i.e. a = 7 = 0.8 , (87) 

which is somewhat smaller than the infrared exponent found above in the angular approximation 

m. 
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The above given infrared analysis is independent of the so far unfixed renormalization constants 
C, and (Recall that the coupled set of DSEs fl2U]) . fIFI]) . flS^ do not depend on From 

eq. flMl) we obtain (x(^) = x(^) + 

lim (a;(A;) - = c , c = ^ - • (88) 

If one uses the infrared expressions for u!{k), d{k),x{k) {and f{k) = 1) defined by eq. (1741) . one 
finds that the energy density is minimized for c = (see ref. [20]). Using the representation [1] 

det J(a) = exp j AxA^ , (89) 

which is correct to the order considered in the present paper, c = imphes an infrared hmit of 
the wave functional 

*(A) = const Yl ^{k) , ^(A; ^ 0) = 1 . (90) 

k 

This wave functional describes a stochastic vacuum, where the infrared modes of the gauge 
field are completely unconstrained. 



5 Numerical results 

The coupled DSEs f fT6|) . fl20|) . fl22l) . fl52l) were solved numerically in the whole momentum range 
as described in refs. [3j and [21j. The renormalization constant was fixed by implementing 
the horizon condition d~^{k = 0) = 0. Furthermore, like in the D = 3 + 1 dimensional case in 
order to stay consistently in 1-loop approximation we have solved the equation for the Coulomb 
form factor f{k) by assuming a bare ghost form factor d{k) = 1 in the DSE for f{k). The 
numerical results obtained are presented in figs. 1 and 2. Fig. 1 shows the gluon energy uj{k) 
and the curvature x(^)- Both quantities are infrared divergent and approach each other for 
A; — ^ in agreement with our infrared analysis given in the previous section. Fig. 2 shows 
the ghost and the Coulomb form factor. The ghost form factor is of course infrared divergent, 
since as shown in sect. 4, self-consistent solution of the Dyson-Schwinger equations exist only 
for infrared divergent ghost form factors. The Coulomb form factor is infrared finite and 
approaches asymptotically for k —>■ oo the ghost form factor. Fig. 3 shows the static color 
Coulomb potential defined by [3] 



where Jo(^) is the (ordinary) zero'th order Bessel function. For the infrared behavior obtained 
in the angular approximation P = 1/2, see eq. (17il) . V{k) ~ 1/^3, which leads to a strictly 
linearly rising quark potential at large distances, while for k oo, where d{k oo) = f{k 
oo) = 1 the potential behaves as V{k) ~ 1/^2 and we obtain the familiar Coulomb potential 
V{r) ~ In^/^ in D = 2 + 1 dimensions. The infrared analysis carried out without resorting 
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Figure 1: The gluon energy u!{k) and the curvature x(A;) obtained from the numerical solution 
of the DSEs for ^ = 2.0(/2 and = 0.5g^. 




Figure 2: The ghost form factor d{k) and the Coulomb form factor f{k) for ^ = 2.0g^ and 
x(/i) = 0.5g\ 
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Figure 4: Comparison of the ghost form factor d{k) obtained from the numerical solution of 
the DSEs for ^ = 2.0g^ and xif^) = 0.5g^ with the lattice data obtained in [9]. 
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Figure 5: The gluon energy uj{k) obtained from the numerical solutions of the DSEs for 
^ = 2.0(7^ and = 0.5g'^ and the corresponding lattice data obtained in [9]. 




Figure 6: The Fourier transform of the static potential obtained from the numerical solution 
of the DSEs and the corresponding lattice data obtained in [9] . 
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to the angular approximation yields V{k) ~ Va;^-^5 k —>■ 0. A careful analysis of our numerical 
solutions (obtained without the angular approximation) yields V{k) ~ Va;^-^5 k —>■ 0, which is in 
between the analytical results obtained with and without the angular approximation. In figs. 
4, 5 and 6 we compare our numerical results for the ghost form factor d{k), the gluon energy 
uj{k) and the form factor of the Coulomb potential (P'{k)f{k) with the lattice data. It has been 
shown by solving the Dyson-Schwinger equation in Landau gauge on the torus [12], [13], that 
very large lattices are required to capture the correct infrared behavior of the Greens functions 
of the continuum theory. While these lattice sizes can be reached in 3 dimensions, they are 
out of reach in D = 4. The lattice calculations performed in ref. [9] in 2 + 1 dimensions used 
lattices of the size 64^, which should be sufficient to extract the correct infrared limit of the 
Greens functions. Our numerical results obtained by solving the DSE in Coulomb gauge are 
in quite satisfactory agreement with the lattice data. In particular, the asymptotic ultraviolet 
and infrared behaviors are quite well reproduced. 



6 Summary and Conclusions 

We have performed a variational solution of the Yang-Mills Schrodinger equation in Coulomb 
gauge in D = 2 + 1. The Dyson-Schwinger equations resulting from the minimization of the 
vacuum energy density have been solved analytically in the ultraviolet and in the infrared and 
in addition some rigorous results of their properties have been derived. In particular, we have 
shown that the ghost form factor as well as the gluon energy have to be infrared divergent, 
which is different from the 3 + 1 dimensional case where solutions of the DSEs exist with 
these quantities being infrared finite [21|. The static non-Abelian Coulomb potential resulting 
from our numerical solution of the DSE is almost linearly rising. Our numerical results are 
in satisfactory agreement with the existing lattice data. The lattice calculations performed 
in D = 3 + 1 so far use too small lattices to give reliable results for the continuum limit, in 
particular, on the infrared properties of the various Green's functions [12] , [13] . 
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